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The running of the 6-quark mass from LEP data * 
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Next-to-leading order QCD corrections to three jet heavy quark production in e + e~ collisions, including quark 
mass effects, are presented. The extraction of the b-quark mass from LEP data is considered and the first 
experimental evidence for the running of a quark mass is discussed. 
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1. Introduction 

The question of the origin of the masses of 
quarks and leptons is one of the unresolved puz- 
zles in present high energy physics. To answer 
this question one needs to know precisely their 
value. However, quarks are not free and their 
mass has to be interpreted more like a coupling 
than an inertial parameter and it can run if mea- 
sured at different scales. Moreover, in the stan- 
dard model (SM) all fermion masses come from 
Yukawa couplings and those also run with the en- 
ergy. To test fermion mass models one has to 
run masses extracted at quite different scales to 
the same scale and compare them with the same 
"ruler". This way, for instance, one can check 
that in some unified models the 6-quark mass and 
the r-lepton mass, although different at threshold 
energies they could be equal at the unification 
scale. 

However, the running of fermion masses, al- 
though predicted by quantum field theory, has 
not been tested experimentally until now. The 
reason being that for energies y/Q* much higher 
than the fermion mass, m q , the mass effects bc- 
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come negligible since usually they are suppressed 
by m 2 JQ 2 . 

While this argument is correct for total cross 
sections for production of heavy quarks, it is not 
completely true for quantities that depend on 
other variables. In particular, it is not true for 
jet cross sections which depend on a new vari- 
able, y c (the jet-resolution parameter that defines 
the jet multiplicity) and which introduces a new 
scale in the analysis, E c = \/Q 2 y c - Then, for 
small values of y c there could be contributions 
coming like m 2 /E 2 = {m 2 /Q 2 )/y c which could 
enhance the mass effect considerably. In addi- 
tion mass effects could also be enhanced by loga- 
rithms of the mass. For instance, the ratio of the 
phase space for two massive quarks and a gluon 
to the phase space for three massless particles 
is 1 + 9,{m 2 JQ 2 ) log(ro g /Q). At Q 2 = m| and 
for the bottom quark this gives a 7% effect, for 
m& = 5 GeV and a 3% effect for = 3 GeV. 

The high precision achieved at LEP makes 
heavy quark mass effects relevant. In fact, they 
have to be taken into account in the tests of the 
flavour independence of a s (mz) This in 

turn means that mass effects have already been 
seen. One can reverse the question and ask about 
the possibility of measuring the mass of the bot- 
tom quark, m&, at LEP by assuming the flavour 
universality of the strong interactions. 

In H we showed that mass effects in three-jet 
production at LEP are large enough to be mea- 
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sured. The observable proposed as a means to 
extract the bottom-quark mass from LEP data 
was the ratio 



rt-(y c )/r 6 



rf,(2/ c )/r d 



(i) 



In this equation T^(y c )/T q is the three-jet frac- 
tion and q denotes the quark flavor. In this ratio 
and at the leading order (LO) the quark mass ef- 
fects can be as large as 1% to 6%, depending on 
the value of the mass and on the jet-resolution 
parameter, y c . 

The three jet decay width is defined by the jet- 
clustering algorithms (see e.g. Q). The effect of 
hadronization is, in principle, small and has been 
estimated using the Montecarlo approach |||| . 

Since the measurement of R^ d is done far away 
from the threshold of 6-quark production it can 
be used to test the running of a quark mass as 
predicted by QCD. 

However, as we discussed in [gj, the leading or- 
der calculation does not distinguish among the 
different definitions of the quark mass: perturba- 
tive pole mass, M b , running mass at M&-scale, or 
running mass at m^-scale. As the numerical dif- 
ference is significant when the different definitions 
of masses are used in LO calculations, in order 
to correctly take into account mass effects, it is 
necessary to perform a complete next-to-leading 
order (NLO) calculation of three-jet ratios includ- 
ing quark masses. 

Although, heavy quark production has been 
considered in a large variety of processes and, 
in particular, in Z-boson decays [|7|-]l0| , there are 
very few NLO calculations of heavy quark jet pro- 
duction taking into account complete mass effects 
( In [jllj and [jiJJ this was done for gluon-gluon 
fusion and virtual-photon production of heavy 
quarks). Only very recently NLO calculations 
of heavy quark jet production in e + e~ collisions 
have become available [p~3| — |57t|] - Here we overview 
our calculation. Its final results were presented 
in 0] and have been used by the DELPHI Col- 
laboration to extract the 6-quark mass at the mz 
scale |||(| showing clearly that indeed the 6-quark 
mass runs from /i = m b to /i = mz as predicted 
by QCD. 



2. Jet ratios with heavy quarks at NLO 

The decay width of the Z-boson into three jets 
with a heavy quark can be written as follows 



r b _ mzg a 

L 3j ~ 



2 647r 2 i9v H v{y c ,r b ) + g A H A (y c ,r b )] ,(2) 



where g is the SU(2) gauge coupling constant, cw 
and sw are the cosine and the sine of the weak 
mixing angle, gv = — I + 4/3s^ and g A = 1 
are the vector and axial- vector coupling of the Z- 
boson to the bottom quark and a s is the strong 
coupling constant. Functions Hy^(y c , r b ) con- 
tain all the dependences on y c and the quark 
mass, r b = (Mb/mz) 2 , for the different algo- 
rithms. These functions can be expanded in a s 
as 



H v(A) =A^(y c ) + ^A^ (A Jy c ) 



-n 



D 



(o) 

V(A) 



(3) 



-B 



(i) 

V{A) 



(Vc n) 



Here, A", is the tree-level contribution in the 
massless limit. It is the same function for the 
vector and the axial- vector parts and it is known 
for the different jet-clustering algorithms in an- 
alytic form. The function ^.y^) gives the QCD 
NLO correction for massless quark. This function 
is also known for the different jet-clustering algo- 
rithms f§ |. Note that even in the chiral limit 
these corrections are different for the vector and 
axial- vector parts. However, the difference, which 
is due to the one-loop triangle diagrams is 
rather small. The net effect of these triangle di- 
agrams is smaller than I0~ 4 in R^ d . In the ratio 
i?3 ( , which is similar to the one defined in eq. ([!]) 
but normalized to the sum of three light flavours, 
l=(u,d,s), the triangle anomaly produces a shift 
of +2 • I0~ 3 . This has been taken into account in 
the experimental analysis Functions 
take into account residual mass effects, once the 
leading dependence in r b has been factorized. The 
tree-level contributions, By^ A y were calculated 
numerically in Q for the different algorithms and 
results were presented in the form of simple fits 

§With our choice of the normalization A^(y c ) = A(y c )/2 
and Ay (y c ) = B(y c )/4, where A(y c ) and B(y c ) are de- 
fined in 
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to the numerical results. Finally, the functions 
By, A y contain the NLO corrections depending 
on the quark mass. 

Note that the way we write Hy^) m ec i- © is 
not an expansion for small rt,. We keep the exact 
dependence on rt, in the functions By/A)- Factor- 
ing out rt, makes it easier to analyze the massless 
limit and the dependence of the results on rj, in 
the region of interest. This means that our results 
can also be adapted, by including the photon ex- 
change, to compute the e + e _ -cross section into 
three jets outside the Z-peak at lower energies or 
at higher energies for top quark production. 

At the NLO we have contributions to the three- 
jet cross section from three- and four-parton fi- 
nal states. For diagrams with emission of four 
real quarks, that can mix different flavours, we 
take the convention that the flavour is defined 
by the quark coupled directly to the Z-boson. 
Therefore, events with emission of a heavy quark 
pair radiated of a light qq are classified as light 
events despite the fact that b quarks are present in 
this four-fermion final state. From the theoretical 
point of view this avoids the appearance of large 
logarithms on the quark mass. The same con- 
vention was considered in the experimental anal- 
ysis H therefore allowing for a consistent com- 
parison. 

One-loop three-parton amplitudes are both 
infrared (IR) and ultraviolet (UV) divergent. 
Therefore, some regularization procedure is 
needed. We use dimensional regularization for 
both IR and UV divergences, because it preserves 
the QCD Ward identities. The three-parton tran- 
sition amplitudes can be expressed in terms of 
a few scalar one-loop integrals Jig] , The result 
contains poles in e = (4 — -D)/2, where D is the 
number of space-time dimensions. Some of the 
poles come from UV divergences and the other 
come from IR divergences. The UV divergences, 
however, are removed after the renormalization 
of the parameters of the QCD lagrangian. After 
that we obtain analytical expressions, which con- 
tain terms proportional to the IR poles and finite 
contributions. 

The four-parton transition amplitudes are also 
IR divergent. These IR divergences cancel the 



corresponding IR poles coming from the virtual 
corrections according to JhJ . Two different meth- 
ods of analytic cancelation of IR singularities 
have been developed: the phase space slicing 
method |2^] and the subtraction method p|, |2l| , p2f . 
We follow the first approach. The four-parton 
transition probabilities for Z — > bbgg(bbqq) are 
split in two parts. The first part contains the 
terms which are divergent when one gluon is 
soft or two gluons (or light quarks) are collinear. 
These terms are integrated analytically in arbi- 
trary D dimensions in the soft and collinear re- 
gions of phase space. This way we obtain the 
IR singular contributions of four partons in the 
three-jet region and show that they are canceled 
exactly by the tree-parton contribution. The 
second part, corresponding to the radiation of 
hard gluons, gives rise to finite contributions and 
can be calculated in D = 4 dimensions. The 
three-jet Z-width is obtained by integrating both, 
renormalized three-parton contribution and four- 
parton transition probabilities, in the three-jet 
phase-space region defined by the different jet- 
clustering algorithms. This quantity is infrared 
finite and well defined. 

Following Ellis, Ross and Terrano fl| (ERT) 
we have classified both, three-parton and four- 
parton transition probabilities, according to their 
color factors. It is clear that the cancelation of 
IR divergences between three-parton and four- 
parton processes can only occur inside groups of 
diagrams with the same color factor. The cance- 
lation of IR divergences can be seen more clearly 
by representing the different amplitudes as the 
different cuts one can perform in the three-loop 
bubble diagrams contributing to the Z-boson self- 
energy. After summing up the three-parton and 
four-parton contributions to the three-jet decay 
width of the Z-boson we obtain the functions 
Hv(A) m eq< @ a t order a s . Since a large part 
of the calculation has been done numerically, it 
is important to have some checks of it. We have 
performed the following tests: i) We have checked 
our four-parton probabilities in the massless limit 
against the amplitudes presented in ref. |2l[] . The 
three-parton amplitudes for massive quarks can- 
not be compared directly with the corresponding 
massless result as they have different structure of 
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IR singularities, ii) The four-parton transition 
amplitudes have also been checked in the case 
of massive quarks by comparing their contribu- 
tion to four-jet processes with the known results 
10 1 . iii) To check the performance of the nu- 
merical programs we have applied our method to 
the massless amplitudes of ERT and obtained the 
known results for the functions iv) We have 

checked, independently for each of the groups of 
diagrams with different color factors, that the fi- 
nal result obtained with massive quarks reduces 
to the massless result in the limit of very small 
masses. 

The last test is the main check of our calcula- 



tion. We have calculated the functions H 



V(A) 



for 



several values of r b , in the range Mb ^1 — 5 GeV, 
and then we have extrapolated the results for 
rb — > 0. In that limit we reproduce the values 
for the function in the different algorithms 
considered and the different groups of diagrams. 
This check is not trivial at all since the structure 
of IR divergences for massive quarks is quite dif- 
ferent from the case of massless quarks: for mas- 
sive quarks collinear divergences are regulated by 
the quark mass, and therefore some of the poles 
in e that appear in the massless case are softened 
by logr & . 

Combining eq. (|l|), eq. (§) and eq. (||) and using 



the known expression for r b m we write R 
the following expansion in a 



as 



R b 3 d = l + r b [b + — h 

7T 



(4) 



where the functions bo and b\ are an average of 
the vector and axial- vector parts, weighted by 

c v = 9vl (9v+9a) and C A = 9a/(9v+9a) respec- 
tively. They can be written in terms of the dif- 
ferent functions introduced before, eq. (||), |§,[l3) 
and also depend on y c and r&. 

It is important to note that because the par- 
ticular normalization we have used in the defini- 
tion of i?3 d , which is manifested in the final de- 
pendence on cy and ca, most of the electroweak 
corrections cancel. Those are about 1% [^3| in 
total rates, while in R\ d are below 0.05%. There- 
fore, for our estimates it is enough to consider 
tree-level values of gv and gA- The same argu- 
ment applies for the passage from decay widths 



to cross sections. Contributions from photon ex- 
change are small at LEP and can be absorbed in 
a redefinition of g v and g\ |J . They will add a 
small correction to our observable. 

Although intermediate calculations have been 
performed using the pole mass, we can also re- 
express our results in terms of the running quark 
mass by using the known perturbative expression 
Ml = mg(0)[l + 2a s (/x)/7r (4/3 - log(m 2 /,x 2 ))] . 
The connection between pole and running masses 
is known up to order a 2 , however consistency of 
our pure perturbative calculation requires we use 
only the expression above. We obtain 



R 



bd 



i+nijj) ( b 



b% - 2b Q log ■ 



,(5) 



where rt(fi) — ml(p)/m 2 z and 

h =6 1 +6 [8/3-21og(r b )] . (6) 

ff,(/i) can be expressed in terms of the running 
mass of the 6-quark at [i = mz by using the 
renormalization group. At the order we are 
working fb(p) = f b (m z ) (a s (m z )/a s (^)y 4:1 °^ 
with a s (fJ.) — a s (mz)/(l + a s (mz)(3ot) and t = 
log(/x 2 /m|)/(4^), /3o = 11 - 2JV//3, N f = 5 and 
7o = 2. 

At the perturbative level eq. (Q) and eq. (||) 
are equivalent. However, they neglect different 
higher order terms and lead to different answers. 
Since the experiment is performed at high en- 
ergies (the relevant scales are mz and mz^/Vc) 
one would think that the expression in terms of 
the running mass is more appropriate because the 
running mass is a true short distance parameter, 
while the pole mass contains in it all the com- 
plicated physics at scales fi ~ M&. Moreover, by 
using the expression in terms of the running mass 
we can vary the scale in order to estimate the er- 
ror due to the neglect of higher order corrections. 
In any case, if one would use in eq. (||) scales as 
low as fi — 5 GeV, one would get something closer 
to the pole mass result. 

Although we have studied the observable 
eq. ([[]) for the four jet-clustering algorithms dis- 
cussed in |||I|,|l3| , in the following we concentrate 
only on the DURHAM algorithm , which gives 
smaller radiative corrections and was the one used 
by the DELPHI collaboration in its analysis. 
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The function bo gives the mass corrections at 
the leading order. As shown in j|] it depends very 
mildly on the quark mass in the region of interest 
( Mb ~ 3 — 5 GeV). Therefore it is appropriate to 
present our results for bo as a fit in only y c : bo = 
En=o k o n} lo g n 2/c- For the DURHAM algorithm, 
in the range 0.01 < y c < 0.10 and 3 GeV < M b < 



5 GeV, using s 



-10.521 , ft, 



(i) 



2 

W 



0.2315, we obtain k, 



(0) 



-4.4352 , ft, 



(2) 



-1.6629 



The function b\ (Tj] gives the NLO massive 
corrections to It is important to note that 

bi contains significant logarithmic corrections de- 
pending on the quark mass. We take them 



into account by using the form 6i = k\ 



(0) 



log(y c ) + km' log(r&) in the fit. The coeffi- 
cients we obtain, for the DURHAM scheme and 
ranges for y c and r& mentioned above, are: k^ = 



.(0) 



297.92 , k[ 1} = 59.358 , k%' = 46.238 



(0) 




Figure 1. NLO results for R b 3 d (DURHAM) for 
(i = mz (dashed), \i = 30 GeV (dashed-dotted) 
and fi — 10 GeV (dotted) for fhb(mz) = 3 GeV 
and a s (mz) = 0.118. For comparison we also 
plot the LO results for M& = 5 GeV (lower solid 
line) and fh^imz) = 3 GeV (upper solid line) 



In fig. [l] we present R\f for fi = mz (dashed) , 
H = 30 GeV (dashed-dotted) and fi = 10 GeV 
(dotted) for fhb(mz) — 3 GeV and a s (mz) = 
0.118. For comparison we also present the LO re- 
sults for the quark mass equal to 5 GeV (lower 
solid line) which is, roughly, the value of the pole 
mass obtained at low energies and 3 GeV (upper 
solid line) which is, roughly, the value one obtains 
for the running mass at the mz scale by using 
the renormalization group. Note that choosing a 
low value for ^ makes the result closer to the LO 
result written in terms of the pole mass, while 
choosing a large /i makes the result approach to 
the LO result written in terms of the running 
mass at the mz scale. 

3. fhb{mz) from LEP data 

If i?3 d is measured to good accuracy one could 
use eq. (||) and the relationship between fhb(fi) 
and fhb(mz) to extract m b (mz)- However, the 
extracted result will depend on the scale \i. For il- 
lustration, in fig. H we show the value for fhb(mz) 
which one would obtain from R b ^ exp = 0.96 (y c 
0.021 as a function of the scale fi. Although one 
would naturally think that the scale /i has to be 
taken of the order ~ mz /3 if the energy is equally 
distributed among the three jets, strictly speaking 
the precise value of \x is undefined. The spread of 
the result due to the variation of the scale in an 
appropriate range gives an estimate of the uncer- 
tainty due to higher order corrections. From fig. |^ 
we see that if we vary /x in the range mz/10 — mz 
the uncertainty in the determination of rhb(mz) 
would be of about 0.20 GeV for the DURHAM al- 
gorithm. In the same range and for the JADE and 
EM algorithms the obtained error would be |13| 
of about 0.25 GeV, while for the E algorithm we 
would obtain an error bigger than 0.50 GeV. 

The DELPHI collaboration has used a slightly 
different approach to extract the value of rhb(mz) 
from (6) . Using eq. (Q) written in terms of the 
pole mass, Mb, and, exploiting the perturbative 
relation between the pole mass and the running 
mass, they obtained the value of fhb(mz) from 
Mb extracted from the experimentally measured 
R\} . The difference between fhb(mz) obtained 
in this way and the running mass obtained di- 
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(i/m 2 



Figure 2. Extracted value of mb(mz) if 

i?3ea;p = 0.96 as a function of the scale /i. We 
take a s (m z ) = 0.118 (solid) and Aa s = 0.003 
(dashed). 



rectly using eq. (|J), which is due to the different 
treatments of the higher order terms, was con- 
sidered as a theoretical uncertainty. Accounting 
in addition for the uncertainty due to the vari- 
ation of the scale /i they get a more conserva- 
tive estimate for the theoretical error, 0.27 GeV , 
for the DURHAM algorithm. We would like to 
note, however, that the central value for rhb(rnz) 
reported in || is fully compatible with the one 
obtained by using only eq. (^|). The compari- 
son of the two methods gives a check of con- 
sistency. Furthermore, the mass of the bottom- 
quark, fhb{mz) 1 measured from the three-jet de- 
cay of the Z-boson || is also fully compatible with 
the value obtained from low energy determina- 
tions p5| after using the renormalization group. 
This provides, for the first time, a nice check of 
the running of a quark mass. 
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